Circulation— Strain Sum Rule in Stochastic Magnetohydrodynamics 
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We study probability density functions (pdfs) of the circulation of velocity and magnetic fields in 
magnetohydrodynamics, computed for a circular contour within inertial range scales. The analysis is 
based on the instanton method as adapted to the Martin-Siggia-Rose field theory formalism. While 
in the viscous limit the expected gaussian behaviour of fluctuations is indeed verified, the case of 
vanishing viscosity is not suitable of a direct saddle-point treatment. To study the latter limit, we 
take into account fluctuations around quasi-static background fields, which allows us to derive a 
sum rule relating pdfs of the circulation observables and the rate of the strain tensor. A simple 
inspection of the sum rule definition leads straightforwardly to the algebraic decay piV) ~ l/F^ at 
the circulation pdf tails. 

PACS: 47.27.Gs, ll.lO.-z 
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I. INTRODUCTION 

Probability distribution functions (pdfs) of quantities 
which "measure" strong fluctuations in a turbulent fluid, 
like velocity differences in Burgers model or the cir- 
culation in three-dimensional incompressible flows (^-|^] 
are a very promising tool in the study of the intermit- 
tency phenomenon. While velocity differences are related 
to the existence of shock waves in the former case, the 
choice of circulation as an observable worth to be investi- 
gated is motivated by the appealing picture of turbulence 
in terms of the complex evolution of vorticity filaments, 
suggested for the first time through direct numerical sim- 
ulations of the Navier-Stokes equations 

An early theoretical analysis of the circulation statis- 
tics by Migdal ||], suggested, with the help of loop equa- 
tion methods, and on the grounds of the central limit 
theorem, that the pdf of the circulation, evaluated for 
a closed contour contained in the inertial range, should 
display gaussian tails. This was not confirmed in direct 
numerical simulations, where pdfs were seen to exhibit 
a stretched exponential decay ^. A further study of 
the problem by Takakura and one of us |^], through 
the instanton approach was able to reproduce ob- 
served results as well as to establish predictions concern- 
ing the role of parity breaking external conditions, which 
presently stand open for experimental verification. 

Our aim here is to address similar questions in the 
problem of turbulent magnetohydrodynamics, motivated 
by the fact that a deeper understanding of the subject 
has been in order for some time, in face of the impressive 
amount of data recorded in astrophysical observations, as 
sunspots and the dynamo effect More specifically, 

we will study the statistics of circulation of velocity and 
magnetic fields in the laminar and turbulent regimes, by 
means of the Martin-Siggia-Rose formalism (MSR) [O. 



The basic technical framework is provided by the saddle- 
point method, where instantons and fluctuations around 
them are assumed to yield, respectively, leading and sub- 
leading contributions to the behaviour of pdfs tails. 

This paper is organized as follows. In sec. II, the basic 
elements of the formalism are set. We define the stochas- 
tic fiuid equations in terms of the Elsasser variables (lin- 
ear combinations of velocity and magnetic fields) from 
which the MSR action and the associated saddle-point 
equations are obtained. In Sec. Ill, as a concrete intro- 
duction to the instanton technique, we discuss the sim- 
pler problem of fluctuations in the viscous limit of the 
magnetohydrodynamical equations, proving circulation 
is governed in this case by exact gaussian statistics, as 
expected on physical grounds. In sec. IV we consider the 
limit of vanishing viscosity. It is known that a direct ap- 
plication of the saddle-point method is fated to fail here, 
as implied from simple dimensional analysis |@,^. We 
circumvent this problem through the strategy devised in 
ref. , which consists of "breaking" the path-integration 
measure into fast and slow degrees of freedom, related, 
respectively, to the vorticity and the rate of the strain 
tensor. The characterization of fast and slow variables 
is the basic physical assumption of our method, as mo- 
tivated by the numerical experiments |^ performed in 
the case of pure turbulence. The outcome of our com- 
putations will be just a sum rule relating the statistics 
of circulation observables and the rate of the strain ten- 
sor. We find then asymptotic expressions for the tails of 
the circulation pdfs, given by the algebraic decay l/F^ 
(the same for velocity and magnetic circulation) , a result 
that indicates the existence of strong intcrmittency ef- 
fects in turbulent magnetohydrodynamics. In sec. V, we 
summarize our findings and discuss directions of further 
research. 
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II. MSR FORMALISM AND THE 
SADDLE-POINT EQUATIONS 

Consider a three-dimensional fluid described by veloc- 
ity and magnetic fields, Va{x, t) and haix, t), respectively, 
which is stirred by large scale stochastic forces. It is pos- 
sible to choose physical units so that velocity and mag- 
netic fields have the same dimensions, and magnetohy- 
drodynamical equations are written in terms of Elsasser 
variables — {va ± ba) @, as 



9tzt + zjdpz^ 

daZ^ = , 



(1) 



where v± — {v±vi,)/2 are the viscosity parameters, which 
account for dissipativc effects, and = {fa ± go) are 
gaussian random external forces, defined through the cor- 
relators 

< >=< >= , 

<ft{s,t)!^{x',t') >= 



25{t-t')[D^^l{x-x')±D%{x-x')] 
5{t~t')Dtp{S~x') , 



(2) 



with 



D'tlix - x')^< Ja{x,t)f0{x' ,t') > 
=< ga{x,t)gp{x',t') > 
= 7?iexp(-|f-ff7LjP) Sap 



D^^l{x-x')=< Ux,t)gp{x',t') > 

= 7?2exp(-|f-ff7L^p) Sap 



(3) 



Above, Li and L2 are the length scales where energy 
pumping mechanisms are assumed to occur, and p > 
parametrizes the spatial decay of the force-force correla- 
tion functions. We note that an extension of the subse- 
quent analysis to alternative definitions of the stochastic 
stirring terms yields no further computational difficulties. 

It is important to determine functionals of the velocity 
and magnetic fields that may carry useful information on 
intermittent behaviour, to be revealed from the tails of 
probability distribution functions. A reasonable choice, 
as discussed in the introduction, is the circulation, de- 
fined in terms of Elsasser variables as 



z^ ■ dx 



(4) 



The above integral is evaluated at time t = 0, assuming 
the fluid evolution started at t = —00. The integration 



contour c is taken to be the circumference 



with z — 0, oriented in the counterclockwise direction, 
i? is a length contained in the inertial range, that is, 
77 ^ R <C Li,L2, where 77 gives the microscopic length 



scale, associated to dominant viscous effects. A further 
physical motivation to study the statistics of (^) comes 
from phenomena usually observed in the astrophysical 
context: due to Ampere-Maxwell's law, the formation of 
magnetic tubes is related to strong fluctuations of veloc- 
ity circulation. On the other hand, by the same token, 
any circulation of the magnetic held leads to the existence 
of conduction/displacement currents, and possible parti- 
cle "jets" , which cross a surface bounded by the contour 
c. 

The joint pdf for the circulation variables may be de- 
fined as 



1 



dX- 



p(r+,r-)^(-)2 

J -00 

X exp(iA+r+ + iX-T-)Z{X+,X-) , (5) 
where the characteristic function is given by 



Z(A+,A-) = 



{exp[—iX'^<j)z^-dx-'iX + (j) zT ■ dx\) 



(6) 



In the following computations, we will consider the an- 
alytical mapping — > iX"^ , so that in the large A^ hmit 
the characteristic function will essentially pick up con- 
tributions from large fluctuations of the circulation (it is 
not difficult to prove this is so for pdfs which decay faster 
than any simple exponential). 

The stochastic partial differential equations (Q) may be 
studied as a field theoretical problem. The MSR formal- 
ism ITlll allows us to write the path integral expression 

Z{iX+,iX-)= J Dz^DDz^DP^DQ^exp{~-S) , (7) 

where the MSR action is 

S = -i J d^xdt[z^{dtz^ + z^dpz^ — v+d^ z^v^c)^ z^ 
+9aP+) + z~{dtz~ + z^dpz^ - v+d'^z^ - v-d'^z^ 

+ daP-) + Q+daZ+ + Q + daZ+] 

-fi f dtd^xd^x'z+{x,t)D+^{x - x')z+{x',t) 



1 



H — / dtd'^xd'^x'z^ {x, t)D^^{x - ^')-^/3 (^'i 



2 _ 

-A+r+ - A^r 



(8) 



A "telegraphic" proof of the above formulation is based 
on the fact that 



1 



exp[ 



dtd-^xd-^xz^ix, t)D^ij{x - x)z'^{x , t)] 



(exp[i / dtd^xz^{x,t)f^{x,t)])f± , 



(9) 



where ((...)}/ stands for the average taken in the ensem- 
ble generated from different realizations of the stochastic 
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forces. Integration over the "response fields" 



yields, 



thus, the fluid equations (0). The role of the and 
is just to assure that daZ^ = daZ^ — 0. In the 
continuum time formulation, it is necessary to introduce 
a jacobian term in (^. However, if time is regarded as 
a discrete variable (as we implicitely may do) there is 
no need for a jacobian (it is unity), usually related to 
anticommuting Grassmann fields |13|| . 

The MSR path-integral formulation may be used to 
recover the Wyld perturbative diagrammatic expansion 
iQ of velocity correlation functions. Since the diagram- 
matic expansion is developed in powers of the convection 
terms, the perturbative MSR- Wyld approach has been 
severely criticized along the years for not taking into ac- 
count singular configurations of the velocity field, which 
are of fundamental importance in turbulence. Neverthe- 
less, one advantage of the MSR formalism is to address 
non-perturbative issues, from the knowledge of specific 
configurations of the flow that give relevant contributions 
to the path integral expression for Z{iX'^ ,iX^). This is 
precisely the task of the saddle-point method, meaningful 
in the limit of large or A~. 

From (^ we get the saddle-point equations 

daZ^^ , 



= , 

+zjdpzt 



dazt 

dtz^ 



dtzt-zjdpz^ + z^dpzj + zjdazj + v+d'^z^ 



i I (f'x D^p{\x ~ x\)z^{x ,t) 
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with 



5T^ 



^ e^p^^Sir^ ^ R)6{z)S{t) 



(10) 



(11) 



where rj_ = {x^ + y'^y^^. 

A major source of difficulty here is a hampering "no- 
go" result that holds in the limit of vanishing viscos- 
ity the saddle-point action computed from solutions of 
( [lO| ) will necessarily depend on A^ in a way incompatible 
with parity symmetry. In fact, a simple check shows that 
saddle-point equations are invariant under a set of scaling 
transformations [^j^, which imply that the saddle-point 
action has the general form 5'^°-' ~ A'^/^, if one takes 
A^ = A and ^ 0. This dependence on A is exactly 
the one found in Burgers turbulence for the statistics of 
velocity differences which we do not expect to re- 

produce, even qualitatively, the parity symmetric pdfs of 
circulation in three-dimensions. In order to find physi- 
cally meaningful results, a solution of this problem was 
advanced in ref. 0, resorting on an alternative formu- 
lation of the MSR path-integral, where the rate of the 
strain tensor is used to parametrize an infinite family of 
saddle-point configurations. We will come back to this 
point in sec. IV. 



III. VISCOUS LIMIT OF CIRCULATION PDFS 

It is interesting to study the viscous limit, where the 
convection terms are neglected in the fluid equations (in 
this case the above no-go theorem does not apply). Let 
us set p = 1 in (||). As a result, we get an instructive 
example where the circulation pdf may be exactly found. 
The saddle-point equations ( |lO| ) are now replaced by 

dtzt+i^+d'^zt - v-d^z^ = 

^-ij d^x'D^f,{\x~x'\)z^{x',t) , 



^^A^ 



SZa 



(12) 



Taking eqs. (g) and (II4), it is possible, with the help 
of some partial integrations, to recast the saddle-point 
action in the simpler form 



S{X, Af,) = —A (p V ■ dx — Xh (p b ■ dx , 



(13) 



where we used A^ = (A ± Af,)/2. It is also convenient 
to deflne Va and ba through z^ — {va ± ba)/2. All we 
need to do, therefore, is to flnd solutions for the velocity 
and magnetic fields at 2: = and t = 0, which we denote 
by Va(x±,0) and ba(xj_, 0), respectively, and to replace 
them in equation ([L3[). 

We may rewrite equations (O) as. 



{dt - vd'^w 
{dt - vbd'')b, 

{dt + iybd^)k 



-I j d''x'D^^l{\x-x'\)vp{x',t) 
-I f d^x-'D^^l{\x-x'\)bp{x',t) 



-tXe3p^^6{r^ - R)S{z)S{t) , 
-i J d^x-'Dl^l{\x-x'\)bp{x',t) 
-I j d'x'D^;^l{\x-x'\)vp{x\t) 

X R 

-iXbCspa — S{r± - R)S{z)S{t) , 



(14) 



Applying {dt + vd^) on the first equation listed in ( p^ ) 
and using also the equations for Va and ba, we will have 

[d'^t - y''{d''f]va{x,t) = ~Fa{x,t) , (15) 

where 

Fa{x,t) = 

= -x( d^x'Di'lQx - x'\)e3^p^S{r'^ - R)d{z')S{t) 



dt + i^bd^' 



)Xb I d'x'D^^l{\x-x'\)e,,^::^S{r'^-R) 



X/ '\xu\ DiX2ttR^ x\ 
xd{z )d{t) ~ e3i3aXi3exip{-j^) 



Ll 



+ ( 



dt + vd^ ,D2Xb2TrR^ 



dt + vbd"^ 



(16) 
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In Fourier space, equation (|16|) becomes 
We obtain, thus, 

1 3 /■ 3-- exp{ik ■ X + iujt) 



Fi'\k) 



+ (- 



iuj — Vfjk- 



■)FL'\k) 



with 



Fa '{k) = -«e3/3afc/3 ^ exp( —) , 



Fi^\k) = -ie^pakp 



i exp(-^) 



(17) 



(18) 



(19) 



Since we are interested to know Va{x±, 0), it foUows, from 



p(r,rb) = ( 



with 



exp 



4C 



4C 



4C 



C = V1V2 - f > 

This condition may be stated as 



V Di 



(27) 



(28) 



(29) 



It follows from (0) and that the above inequality is 
always satisfied (in fact, Di > D2 while the Ihs of ( [29| ) 
is > 2). Considering, for instance, D2 = 0, which gives 
773 = 0, we find the factorized form of the circulation pdf. 



p(r,rfc) = ( 



1 



exp 



4?7i 



4?72 



(30) 



Va{x±,0) = / d^k exp{ik± ■ x±) 



1 

X — 

k^ 



Fi'\k) , Fi^\k) 



Substituting (|l|) in (|o|), we find 



Va{x±,Q) = 



2u {u + Vb) 



Similarly, 



ba{x^,Q) 



ttR 



2 r 



2zyb + i^b) 



(20) 



(21) 



(22) 



Thus, from (|l|), (|l]) and (g), we obtain the saddle- 
point action 



S'(A, Afc) = -A^7?i - A^772 - AAf,773 . 



where 



3^ 



m 



, '72 



, V3 = 



3(i^ + i^b) 



(23) 



(24) 



Performing now A —>■ ~iX and Ab — > — iA;,, to restore the 
original definition of these parameters, we get 

Z(A, Ab) cx exp(-A^7;i - Ab??2 - XXbris) . (25) 

Defining 

(26) 



T = j) V ■ dx , Tb = j) b ■ dx , 
we obtain, from (a), the gaussian pdf 



IV. INVISCID THEORY AND THE 
CIRCULATION-STRAIN SUM RULE 

Let us study now the inviscid theory, where i^, t'b — >■ 0, a 
limit associated to the fully developed turbulent regime. 
Since the radius R of the contour c is much smaller than 
Li and L2, the scales where energy is injected into the 
system, we are allowed to consider linear expressions for 
the fields z^, 



(31) 



with the tensor of velocity derivatives satisfying to 
X^a ^aa = Oj ^lue to the incompressibility constraint and 
the absence of magnetic monopoles. To work in the order 
of approximation given by (|3l|) we consider a quadratic 
form for the pressure fields ^ 



(32) 



so that the gradient terms daP^ exactly cancel in ( [lO| ) 
any symmetric tensor acting on the spatial coordinates, 
which would appear in the linear approximation. We 
are left this way with equations which describe the time 
evolution of the antisymmetric part of a^p (related to 
the vorticity) 



dt'^'^f'' 



f d^xd^^p^^^{\x\)z^{x,t) 



(33) 



where we have defined 



d[^P%^{\x\) = \ (d^D^^{\x\) - dpDt^{\x\)) . (34) 
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Before proceeding, it is necessary to discuss the alter- 
native definition of the saddle-point method that we will 
employ. It has been suggested through numerical experi- 
ments of pure turbulence ||] that the rate of strain tensor, 
a" , does not fluctuate so strongly as the vorticity. As it 
was shown in ref. 0, one may take advantage of this 
physical observation to overcome the no-go result com- 
mented in sec. II. We will assume that the characteriza- 
tion of symmetric and antisymmetric degrees of freedom 
as slow and fast fluctuating variables, respectively, holds 
also in the magnetohydrodynamical realm, where similar 
vorticity filamentary structures are as well observed. Of 
course, experiments will decide ultimately if this assump- 
tion is correct or not, but we provisionally regard it as a 
working hypothesis that allows us to establish meaning- 
ful predictions. Following ref. we may define the MSR 
path-integral for the characteristic functional as 

Z{iX+,i\-)^ j Da^' J Dz^Dz^ 

xDP^DQ^DQ^exp{~S) , (35) 
with a^^ = iT^''{x,y,t) and 



S - - I dxdydtQ'^^{x,y,t) ■ [daz^l^^o + dpz+l^^o 



I 
2 

-2cr+J(a;,y,i)] - ^ J dxdydtQ~p{x,y,t) ■ [daZ'^\^^o 
+a/3Z-U=o-2a-^(x,y,t)] . (36) 
The functional form ( ^5|) is obtained by putting 

1 = / Da^'DQ^ exp(5 - S) , (37) 



in the integrand of (|^) . The order of integrations is inter- 
changed, to write the integral over a^p as the last one to 
be performed. These mathematical steps are motivated 
by the fact that the rate of the strain tensor plays the 
role of a quasi-static background where vorticity fluctu- 
ations take place. The central idea is then to apply the 
saddle-point method to the action S, treating cr^^ as ex- 
ternal fixed fields. The saddle-point equations for z^, in 
(p^), are replaced now by 

dtz^ - zjdazj + zjdpz^ + z^dpz"^ + v+d'^z'^ 



+v^d^z^ + 9„Q± + dp{5{z)Qf,J + iX- 



--^11 

5Za 



0, (38) 



There are also two additional equations, associated to 
variations of the fields Q^^, 

9a2:|U=o + dpz^^a - 2(7±|(a;,y,t) = . (39) 

In order to seek for solutions of the saddle-point eqs. 
(p^), note they are invariant under rotations around the 
z axis. The most general form of an axisymmetric tensor 
of velocity derivatives is given by 



a±(0 



a±(t) uj^it) 
-w±(i) a±(t) 
-2a±(i) 



(40) 



From now on, we will substitute cr^^ appearing in the 
above relations by the axisymmetric expression (|4C 
This approximation amounts to the replacement 



Da^{t) , 



(41) 



in the path- integral (p5|). If one were able to integrate 
exactly the expression for Z(A^), keeping a*(t) fixed, the 
circulation pdf could be written as 



p(r±) = j Da^-p[a^W^\a^] 



(42) 



where p\a^] is the probability density functional to get 
the axisymmetric rate of the strain tensors defined by 
a*(t), and p[r*|a^] is the conditional pdf to get F^, 
in the background a^{t). A natural question is how to 
reproduce this kind of relation through the instanton ap- 
proach. Is the saddle-point action leading to an approx- 
imation for p[a^], p[F^|a*], or both?. We do not know 
how to answer it a priori. We take a pragmatical point 
of view, where an answer is found only after concrete 
computations are performed. 

Using ( prj ) , ( ^ ) , ( |40| ) , and taking the limit of vanishing 
viscosity, we may write (^8|) as 



^i\^ezo.p^J{r^-R)5{z)8{t) . 



(43) 



We have now a closed set of coupled equations given by 
(B3) and \\ 



A crucial observation is that the viscosity 
terms in (|38) have the opposite sign, if compared to the 
ones appearing in the usual fluid equations. We have 
to impose, therefore, in order to avoid an unbounded 

= 0, for t > 0. In 



growing of the fields (x, f) , that 



this way, (MS) leads to the boundary condition 



z±(af, 0-) = a±e3/5a— <5(r± - i?)(5(z) 

r I 



(44) 



Also, we require that z^{x,t) — > as i ^ — oo. The 
equation for z^ {x, t) may be solved through the ansatz 



z^ix,t) = e,(3^xpSiz)J2ci{t)rl-'6(^\r^ 



R) , (45) 



where (5("H^-L - R) ^ d''S{rj_ ~ R)/drl. The problem of 
finding z^ is mapped into the computation of c^{t). The 
boundary condition (44) reads now 



c±(o-) = a± , 

ct{0~) = 0, for n > . 
We obtain, substituting (E^) in (E^), 



(46) 
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-Cn = [A~{2 + n)B]Cn - BC„-i 



(47) 



with C_i = and 

Ait) 



a-(t) 
a+{t) 



Bit) 



a- it) 
a+it) 



Cnit) 

Furthermore, we get Q — 0, and 
'^ir±,t) 



(48) 



(49) 

An exact solution for C„(i) may be found in the case 
where [A, -B] = 0. Matrices A and B commute only if 
a~^it) — it) = ait), which means that 

dabfl + dfjba = , 

daVp + dfsva = 2ait) ■ iS^p - SJsq^s^) . (50) 

We are concerned, thus, with fluctuations of the circu- 
lation variables in the presence of a vanishing rate of 
the strain tensor of the magnetic field. Even though the 
need of computational feasibility compels us to the study 
of the conditional pdf of circulation variables in flow con- 
figurations where (^0|) is verified, this restriction is not 
unfortunate by no means: we will be able, below, to find 
well-defined predictions related to intermittent fluctua- 
tions of the velocity or magnetic circulation. 
The exact solution of (^7|) is given by 

r< n\ * ~ r dt'{2B-A) f -f*dt'B lA^T rt^i\ 
Cnit) ^ —e Jo Me Jo -11 A, (51) 



where 



(52) 



Taking into account (^), the inflnite series ( ^ ) may be 
exactly summed up, to yield, in terms of Va and ba, 

Vaix,t)^ iXe3fja—Sirj_ - Re-fo ^)Siz) , 



bnix,t)= i\he^Br, — e Jo 



x,5(ri-i?e/o'^*''^(*'))5(z) . 



(53) 



To find the saddle-point action S*^"-*, it is necessary to get 



-it). Using ([40|), we write eqs. ( |33| ) as 



d + ,4 
— a; + aiu!^ 



= -ij d^xd[^p^^^i\x\)ztiS,t) . 



(54) 



Substituting the solutions for (z, t) in (| 
for t < 0, 



d 



we obtain, 



— + aiuj^ 



where 



r;± , ;± -2 r*<it'a(t')n 2pfdt'a(t') 
Vl + '2 S " ^ '\e Jo 



(55) 



(56) 



The idea now is to consider, in the context of a gradient 
expansion, the effects of time independent configurations 
ait) = a. Therefore, we replace the path-integration over 
arbitrary fields ait) by an ordinary integration over a. 
Furthermore, it turns out that w*(i) — > as t — > — oo 
only if a > and p > 1. Provided these conditions are 
satisfied, we obtain 



o.±(0 = 



i2p + \)if-q 



+ 



4ap(;, + Ty^^"P[2pat] 

^ exp[2(p - l)at] 



i2p~ l)lf ~ 



Aapip — 1) 

The saddle-point action may be written as 



(57) 



2 



dt 



d^xd^x' 



x[z+(x, t)D+f^ix-x')z+ix,t) 
+z-ix, t)D^pix - x')z^ix\ t)\ 



-(0) • dx 



(58) 



Using ( |53| ) and ( |55| ) , a straightforward computation gives 
a quadratic form in A and A;, for the action S^^K With 
the help of the notation introduced in the analysis of the 
viscous limit, equation (p3), we have now 



Dipfjp) 
Ll'ip+l) 



V2 = 



DiPfjp) 
L?ip-l) 



2D2fip) 



-2p 



where, in terms of gamma functions. 



fip) 



7r2i?2(p+i) 



r(2p+l) 



rip + i)Tip + 2) 



(59) 



(60) 



The requirement of having a positive definite saddle- 
point action S**^"^ (after the analytical mapping A^ 
iX^) leads to /(p) > 0, that is p < 2, and 



4>i 



D2L'/ 



(61) 
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Therefore, our formalism is assumed to work in the in- 
viscid theory for 1 < p < 2, if the above condition is 
also verified (what happens, in particular, for D2 = 0). 
It is interesting to observe that 772 diverges when p — > 1, 
so that the pdf of the magnetic circulation becomes very 
broad in that limit. This fact may be physically inter- 
preted as due to a progressive decoupling of magnetic and 
velocity fields as p approaches unity, in flow realizations 
where the magnetic strain vanishes. The velocity field 
would then behave in the same fashion as in pure hy- 
drodynamical turbulence, while the magnetic circulation 
would diffuse in a strong way, under the action of the 
external stochastic forces. On the other hand, for p large 
enough, higher wavenumbers in the Fourier transform of 
the force-force correlation functions cannot be neglected, 
and the saddle-point method, as based on smooth instan- 
ton solutions, breaks down. 

Taking the definitions for the ?7's in (|59|), the condi- 
tional pdf for the circulation variables, p(r,rf,|a), may 
be readily obtained from the expression previously de- 
fined in (P?!). We may write now, recalling (^), the joint 
circulation pdf as 



p(r,rfc) 



dap{a)p{r,Tb\a) 



(62) 



This is a sum rule which holds for asymptotically large 
values of F or F;,. It is important to keep in mind that 
the pdfs appearing in (|6^) are defined under the condition 
that dabfs + dpha = 0. Also, it is natural to expect that 
the unknown function p{a) has a finite limit as a ^ 0, so 
that we are led to the algebraic decay at the pdf tails 



p(r,Ffc)^ 

Dip 



Dip 



Lfip^l) 



Ll'ip+l) 



2D2 

j2p 



FFft 



-3/2 



(63) 



a result which clearly signalizes the intermittent nature 
of circulation fluctuations in turbulent magnetohydrody- 
namics. In particular, we may get from ( p3|) the pdf tail 
for the magnetic circulation as 



p(Ffc) = j dFp(F,Ffc) ^ iFfcl- 



(64) 



An analogous behaviour follows for the velocity circula- 
tion pdf. Since the magnetohydrodynamic system is most 
of the time around the state where the magnetic strain 
vanishes (that is, the mean magnetic strain is zero), we 
conjecture that ( ^4| ) is a general result, holding beyond 
a specific selection of ensembles in the fully turbulent 
regime. It is worth mentioning that similar computations 
for the case of pure turbulence reveal "less intermittent" 
fluctuations, i.e, the circulation pdf tail p(F) ~ l/|Fp. 

One could be puzzled by the fact that we have found 
algebraic decaying pdfs, exploring the large limit, 
which, as commented before, is related to pdfs with fast 
exponential decay. However - a crucial point in the anal- 
ysis - the instanton method becomes meaningful since 



it has been used to get the conditional pdf p(F,Fb|a), 
which has an exact gaussian shape. We stress that the 
algebraic decay in eq. (|6^) follows from the integration 
over the velocity strain parameter a in eq. (|6^). Also, at 
this point it should be noted that the prediction of the 
stretched exponential tail for the circulation pdf in pure 
turbulence, as previously discussed in ref. 0|, is indeed 
likely to hold at an itermediate range where deviations 
from gaussian statistics appear, and not at the extreme 
asymptotic region. To describe the intermediate range - 
the one so far observed in direct numerical simulations - 
it is necessary to know in some detail the form of p{a), 
the strain pdf. 



V. CONCLUSIONS 

We studied the problem of stochastic magnetohydro- 
dynamics in the limits of large and small viscosity pa- 
rameters, focusing on circulation variables, which are ex- 
pected to yield information on itermittent behaviour in 
the latter situation (turbulent regime) . From the techni- 
cal point of view, the computational framework was given 
by the application of the saddle-point method within the 
Martin-Siggia-Rose path-integral formalism. 

The instanton approach was straightforwardly applied 
to the viscous case, where we found a gaussian form for 
the circulation pdf, as it should be. Regarding the in- 
viscid limit, when turbulence surely comes into play, our 
investigation is particularly devoted to the properties of 
pdf tails. The fundamental physical hypothesis is that for 
both velocity and magnetic fields, the rate of the strain 
tensor behaves as a quasi-static background where faster 
fiuctuations of the antisymmetric part of the tensor of 
field derivatives (related to the vorticity) occur. This mo- 
tivates a definition of the Martin-Siggia-Rose functional 
which explicitcly takes into account fast and slow degrees 
of freedom. The saddle-point method is then applied to 
the action defined in terms of the fast variables. Such a 
strategy gives a way out of the problematic dimensional 
constraints imposed on the form of the saddle-point ac- 
tion in the original path-integral formalism. We were able 
in this way to obtain a sum rule relation for the statis- 
tics of the circulation variables, as well as a quantitative 
prediction for the algebraic decay of pdf tails. 

Some mathematical simplifications were employed in 
the course of analysis, which basically fall into two 
classes: either based on the role of fast and slow vari- 
ables or related to the manipulation of exact solutions of 
the saddle-point equations. Improvements on the latter 
aspect are likely to be the more interesting (and per- 
haps the more difficult to be attained). The relevant 
open question is then how to implement the instanton 
approach when exact saddle-point solutions are not avail- 
able anymore, in order to compute joint pdfs of circula- 
tion variables in turbulent magnetohydrodynamics. 
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